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We analyze a gap equation for the propagator of Dirac quasiparticles and conclude that in 
graphene in a magnetic field, the order parameters connected with the quantum Hall ferromag- 
QO ' netism dynamics and those connected with the magnetic catalysis dynamics necessarily coexist (the 

latter have the form of Dirac masses and correspond to excitonic condensates). This feature of 
graphene could lead to important consequences, in particular, for the existence of gapless edge 
states. Solutions of the gap equation corresponding to recently experimentally discovered novel 
^ \ plateaus in graphene in strong magnetic fields are described. 
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The properties of graphene, a single atomic layer of graphite^, have attracted great interest, especially after the 
experimental discovery- 2 * 3 - and (made independently) theoretical predictions^ of an anomalous quantization in the 
quantum Hall (QH) effect. In this case, the filling factors are v — ±4(|n| + 1/2), where n is the Landau level index. 
For each QH state, a four-fold (spin and sublattice-valley) degeneracy takes place. These properties of the QH effect 
are intimately connected with relativistic like features in the graphene dynamics 

In recent experiment s 10 ' 11 , it has been observed that in a strong enough magnetic field, B > 20 T, the new QH 
plateaus, v = 0, ±1 and ±4, occur, that was attributed to the magnetic field induced splitting of the n = and 
n = ±1 Landau levels (LLs). It is noticeable that while the degeneracy of the lowest LL (LLL), n = 0, is completely 
lifted, only the spin degeneracy of the n = ±1 LL is removed. 

On theoretical side, there are now two leading scenarios for the description of these plateaus. One of them is 
the QH ferromagnetism (QHF)— ' 13 ' 14 ' 15 (the dynamics of a Zeeman spin splitting enhancement considered in RefJ£ 
is intimately connected with the QHF). The second one is the magnetic catalysis (MC) scenario in which excitonic 
condensates (Dirac masses) are spontaneously produce d 17 ' 18 ' 19 ? 20 . For a brief review of these two scenarios, see Ref3I. 

While the QHF scenario is based on the dynamical framework developed for bilayer QH systems 2 ^, the MC scenario 
is based on the phenomenon of an enhancement of the density of states in a strong magnetic field, which catalyzes 
1— — '■ electron-hole pairing (leading to excitonic condensates) in relativistic like systems. The essence of this effect is the 
\ dimensional reduction D — > D — 2 in the electron-hole pairing dynamics and the presence of the LLL with energy 
E = (containing both electron and hole states) in relativistic systems in a magnetic field. This universal phenomenon 
f* — ' was revealed in Ref3^ and was first considered in graphite in Refsj 24 ' 25 . 

On technical side, the difference between these two scenarios is in utilizing different order parameters in breaking 
the spin-sublattice-valley £7(4) symmetry of the noninteracting Hamiltonian of graphene. While the QHF order 
parameters are described by densities of the conserved charges connected with diagonal generators of the non-abelian 
1 subgroup SU(4) C £7(4), the order parameters in the MC scenario are Dirac mass like terms. Note that while the 
latter are bifermion operators which are invariant under 2 + 1 dimensional Lorentz transformations (with the Fermi 
velocity vf — 10 6 m/s playing the role of light velocity), the QHF charge densities are time like components of the 
corresponding conserved currents which transform as vectors under the Lorentz transformations. 

One may think that the QHF and MC order parameters should compete with each other. However, as will be 
shown in this paper, the situation is quite different: These two sets of the order parameters necessarily coexist, which 
implies that they have the same dynamical origin. The physics underlying their coexistence is specific for relativistic 
like dynamics that makes the QH dynamics of the £7 (4) breakdown in graphene to be quite different from that in 
bilayer QH systems^ whose dynamics have no relativistic like features. 

Our approach is based on studying the gap equation for the propagator of Dirac quasiparticles. For the description 
of the dynamics in graphene, we will use the same model as in Refsi^ 4 ^, in which while quasiparticles are confined to 
a 2-dimensional plane, the electromagnetic (Coulomb) interaction between them is three-dimensional in nature. The 
dynamics will be treated in the Hartree-Fock (mean field) approximation, which is conventional and appropriate in 
this cas o 12 ' 13 ' 17 ' 24 ' 25 . Then, at zero temperature and in the clean limit (no impurities), the gap equation takes the 
form: 

G-\x, y) = S-\x, y) + ih 7 °G(x, yh°6(x - y )U c (x - y) 

-ihrftotfGfaxfiS^x - y)t^(O). (1) 

Here x = (xq,x), with xq = t being time coordinate, Uc(x) is the Coulomb potential in a magnetic field, given in 
Eq. (46) in Ref. 25 , £7^(0) is its Fourier transform at k = 0, G^ 1 (x,y) is the full inverse quasiparticle propagator, and 
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S (x,y) is the bare inverse quasiparticle propagator, 

iS-\x, y) = [(ihd t + ^ - fiBBa 3 )^ - v F ^{x - y), (2) 

where /zo is the electron chemical potential, ir = —ihd + eA/c is the canonical momentum, and ^,bB^ <j 3 is the 
Zeeman term [the vector potential A corresponds to the magnetic field B, B = \B\, fiB is the Bohr magneton, and 
the Pauli matrix a 3 acts on spin indices) 2 ^. 

For Dirac matrices 7 , 7, we use the same representation as in Refsi 17 ' 25 (xy-plane is chosen for graphene). Note 
that while the second term on the right hand side of Eq.([T]) describes exchange interactions, the third one is the 
Hartree term describing annihilation interactions. 

The analysis of gap equation |T]) beyond the LLL approximation is a very formidable problem. Because of that, 
we will utilize the following approximation: the Coulomb potential Uc(x) in the gap equation will be replaced by the 
contact interaction Gi n tS 2 (x): 

G-\x,y) = S-\x,y) + ihG mtl °G{x 1 x) 1 8 3 {x - y) - %hG int f ti[ 7 °G(x,x)]S 3 (x - y), (3) 

where Gi n t is a dimensional coupling constant. Such an approximation is common in Quantum Chromodynamics 
(QCD), where long range gluon interactions are replaced by contact (Nambu-Jona-Lasinio) ones. This leads to a good 
description of nonperturbative dynamics in low energy region in QCD (for a review, see for example Ref. 27 ). Because 
of the universality of the MC phenomenon and because the symmetric and kinematic structures of the gap equations 
((T|) and ([3]) are the same, we expect that approximate gap equation © should be at least qualitatively reliable for the 
description of the LLL and first few LLs, say, n = ±1 LL. This in turn implies that in the analysis of this gap equation, 
one should use an ultraviolet cutoff A of the order of the Landau scale L(B) = y/h\eB±\v F /c ~ 300^/£>^(T) [K] (in 
Kelvin), where Bj_ is the component of B orthogonal to the graphene plane measured in Tcsla. The dimensional 
coupling constant Gi nt should be taken then as Gi nt ~ 1/ VeB± (see below 

Because of the Zeeman term, the U(4) symmetry is broken down to the "flavor" symmetry U(2) + x t/(2)_, where 
the subscript ± corresponds to up and down spin states respectively. The generators of the U{2) Sl with s = ±, are 
I <8> P s , — «7 3 <S> P s , 7 5 ® P s , and 7 3 7 5 <g> P s (here I is the 4x4 unit matrix, 7 s = i7°7 1 7 2 7 3 , and P± = (1 ± er 3 )/2 are 
projectors on spin up and spin down states)^. 

Our goal is searching for solutions of Eq. © both with spontaneously broken and unbroken SU(2) S , where SU(2) S 
is the largest non-abelian subgroup of the U(2) s . The Dirac mass term A s 4!P s ip = Ag^^Pgip, where A s is a Dirac 
gap (mass), is assigned to the triplet representation of the SU(2) S , and the generation of such a mass would lead to 
spontaneous flavor SU(2) S symmetry breaking down to the U(l) s with the generator 7 3 7 5 ® P A 17 i 24 i 25 . There is also 
a Dirac mass term of the form A S ^7 3 7 5 P S ^ that is a singlet with respect to SU(2) S , and therefore its generation 
would not break this symmetry. On the other hand, while the triplet mass term is even under time reversal T, the 
singlet mass term is T-odd (for a recent review of the transformation properties of different mass terms in graphene, 
see Refj^2). It is noticeable that consequences of the presence of the mass A in graphite were discussed long ago in 
Ref A 

The analysis of gap equation ^ that we use is closely connected with that in Ref3L and based on the decomposition 
of the quasiparticle propagator over the LL poles with the residues expressed through the generalized Laguerre 
polynomials. A detailed description of the analysis will be presented elsewhere. Here we will describe its main results. 
It was found that, for a fixed spin, the full inverse quasiparticle propagator takes the following general form (compare 
with Eq. (0): 

iG-^x, y) = [{ihdt + Hs+ A^Vh - «f??7 - A s + A s7 3 7 5 ]5 3 (x - y), (4) 

where the parameters fi s , jj, s , A s , and A s are determined from gap equation ©. Note that the chemical potential 
fi± includes the Zeeman energy =fZ, with Z — hbB — 0.67£>(T) [K], and the chemical potential jl s is related to 
the density of the conserved pseudospin charge ip^ 7 3 7 5 P s ip , which is assigned to the triplet representation of the 
SU(2) S . Therefore, while the masses A s and A s are related to the MC order parameters {ip^ 3 ^ 5 P s ip) and (-ipP s ip), the 
chemical potentials /13 = (/i + — fi-)/2 and jl s are related to the conventional QHF ones: the spin density (ip^a 3 ^) and 
the pseudospin density (V'^7 3 7 5 -P s ^), respectively. Note that while the triplet Dirac mass term describes the charge 
density imbalance between the two graphene sublattice o 17 ! 24 , the pseudospin density describes the charge density 
imbalance between the two valley points in the Brillouin zone. 

The dispersion relations for higher LLs (|n| > 1) following from Eq. (T?J) are 



E$ = -fi s + a jl s + sign(n)^2h\neB ± \v 2 F /c + (A s + aA s ) 2 , 



(5) 
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where a = ±1 are connected with eigenvalues of the pseudospin matrix 7 3 7 5 . The case of the LLL is special, and its 
dispersion relation is 

E { s a) = -fi s + a[fi s sign(eS_L) + AJ + A s sign(eS ± ). (6) 

One can see from Eqs. ([5]), © that at a fixed spin, the terms with a are responsible for splitting of LLs. 

In fact, for each value of spin, our analysis revealed the following three types of solutions: a) a singlet solution with 
a nonzero singlet mass A and with no triplet parameters A and fl, b) a triplet solution with nonzero A and /2, and 
with the singlet mass A being zero, and c) a mixed solution with A, A, and fl being nonzero. The latter is realized 
only in higher LLs. In order to find the most stable solution among them, we compare the free energy density f2 of the 
corresponding ground states. In the mean field approximation that we use, f2 takes the following form on solutions of 
the gap equation^ 



QVT = iTr 



LnG*- 1 + - (S^G - 1) 



(7) 



where VT is the space-time volume, the trace, the logarithm, and the product S 1 G are taken in the functional sense, 
andG = diag(G+,G_). 

The process of filling the LLs is described by varying the electron chemical potential /io- We will consider positive 
Ho (dynamics with negative /io is related by electron-hole symmetry and will not be discussed separately). In this 
paper we will mostly consider the LLL dynamics (results for the n = 1 LL will be briefly described at the end of the 
paper). 

For the case when only the LLL is doped, which corresponds to the condition |/i s ± fl s \ <C L(B), we arrive at the 
following results: 

i) A solution with singlet Dirac masses both for spin up and spin down is the most favorable for < /io < 2A + Z, 
where A = G lnt \eBj_\/8Trhc—. It is: 

A± = A± = 0, (jl± = p.± T A, A± = ±M sign(eS ± ) (8) 

with p,± = fio =p Z and M = A/(l — A), A = X^/ttL 2 (B) /2A where the dimensionless coupling constant A is A = 
G int A/ (ATr 3 / 2 h 2 Vp)^: From dispersion relation we find that E + > and £7_ < 0, i.e., the LLL is half filled 
(the energy spectrum in this solution is a independent). Therefore the spin gap AE = E + — E_ corresponds to the 
v = plateau. The value of the gap is AEq = 2M + 2(Z + A). It is instructive to compare AEq with the spin gap 
in Ref^. The latter contains an enhanced Zeeman spin splitting, which corresponds to the second term 2(Z + A) 
in A£? . However, besides this term, there is also the large contribution 2M in AE in the present solution, which 
is connected with a dynamical singlet Dirac mass for quasiparticles. The presence of this mass could have important 
consequences for gapless edge states whose relevance for the physics of the v = plateau was pointed out in Ref^. 
Generalizing the analysis ir>22, we have found that such states exist only when the full Zeeman splitting Z + A is 
larger than the Dirac gap M = A/(l — A). This leads to the constraint Z > A Aj (1 — A). Let us consider the case with 
B = B^. Then, since Z ~ B± and A ~ y/B± (see below), this constraint leads to a lower limit B^^ for the values of 
B± at which gapless edge states exist. On the other hand, since Z depends on total B while A depends only on B±, 
adding a longitudinal B\\ will decrease the lower limit for B±. It would be interesting to check experimentally this 
point. Also, these features could be relevant for the interpretation of the recent experiments^, in which no gapless 
edge states were detected for B = B± < 14 T. We shall return to this issue below. 

ii) A hybrid solution, with a triplet Dirac mass for spin up and a singlet Dirac mass for spin down, is the most 
favorable for 2A + Z < fi < 6A + Z. It is: 

A+ = M, fi+ = Asign(eS ± ), M+ = fi+ - 4A, A + = 0, 



A_ = /i_ = 0, /i_ = /!_ - 3A, A_ = -Msign(eBx). (9) 

As follows from Eq. ©, while E { + 1] > 0, the energies E+ ^ and = E_ ^ are negative. Consequently, the 

LLL is now three-quarter filled and, therefore, the gap AE\ = E+ — E{_ = 2(M + A) corresponds to the v = 1 
plateau. The latter, unlike the v — plateau, is directly related to spontaneous SU(2) + flavor symmetry breaking. 

iii) A solution with equal singlet Dirac masses for spin up and spin down states is the most favorable for /io > 6A+Z. 
It is 



A± = A± = 0, /i± = ji± - 7 A, A± = -Msign(eS_L) 



(10) 
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(compare with Eq. (JSj) ) . It is easy to check from ([6]) that both E+ and £L are negative in this case, i.e., the LLL is 
completely filled. Therefore, this solution corresponds to the v = 2 plateau related to the energy gap AE2 — \2L(B) 
between the LLL and the n = 1 LL. 

This analysis leads us to the picture for the LLL plateaus which qualitatively agrees with that in experiment o 10 ' 11 . 
In particular, taking the dimcnsionlcss coupling A to be a free parameter and choosing cutoff A to be of the order 
of the Landau scale L(B), we arrive at the scaling relations, A ~ \/\eB±\, M ~ y / [eBjJ, and, therefore, AEi = 
2(A+M) ~ ^/\e~B±\ for the gap related to the v = 1 plateau. One can check that the experimental value AEi ~ 100 K 
for B± — 30 T— corresponds to A ~ 0.02. However, because interactions with impurities are ignored in the clean 
limit used in the present model, it would be more reasonable to consider A, say, in interval 0.02 — 0.2. Then, for these 
values of A, we find from the constraint Z > XA/ (1 — A) in the solution i) above that the gapless edge states exist for 
\B±\ > B [ l r \ where 0.01 T < s[ cr) < 200 T. One can see that B'f^ is sensitive to the choice of A. Therefore in order 

to fix the critical value B± r ^ more accurately, one should utilize a more realistic and constrained model. 

As to the n = 1 LL, we found that there are the gaps AE 3 = AE 5 ~ 2A and AE 4 ~ 2(Z + A) corresponding 
to the plateaus v = 3, 5 and v = 4, respectively (the contributions of Dirac masses are suppressed at least by factor 
M 2 /L 2 (B) there). Note that A-E^g and AE4 are essentially smaller than the LLL gaps AE\ and AEo, respectively 
(Ai?3 5 < AE\/2). On the other hand, the experimental data yield AE4 ~ 2Z, and no gaps AE3, AE$ have been 
observe d 10 ! 11 . We believe that a probable explanation of this point is that, unlike Z, the value of the dynamically 
generated parameter A corresponding to the \n\ > 1 LLs will be essentially reduced if a considerable broadening 
of higher LLs in a magnetic field is taken into account^. If so, the gap AE4 will be reduced to 2Z and the gaps 
Ai?3, AE$ will become unobservable. 

Recently, in Refi2£, a large width Ti of 400 K was determined for the n = 1 LL. The plateaus ^ = 3,5 could become 
observable if the gaps AE3 = AE$ ~ 2A calculated in the clean limit are at least of order Ti or larger— . The LLL 
gap AEi ~ 100 K at \B±\ = 30 T corresponds to AE$^ < 50 K. Then, taking a conservative estimate Li = 100 K and 

using A r~j y/\eB±\, we conclude that to observe the v = 3, 5 plateaus, the magnetic fields should be at least as large 
as B ~ 100 T. 

In conclusion, we have shown that the QHF and MC order parameters in graphene are two sides of the same coin 
and they necessarily coexist. This feature could have important dynamical consequences for low energy excitations, 
in particular, for gapless edge states. It would be desirable to extend the present analysis to a more realistic model 
setup, including the genuine Coulomb interactions, LLs impurity scattering rates, and temperature. 
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